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Abstract 

We show, in general, that when a discontinuity of either zeroth-order or first-order takes place in an order 
parameter such as the chiral condensate, discontinuities of the same order emerge in other order parameters 
such as the Polyakov loop. A condition for the coexistence theorem to be valid is clarified. Consequently, 
only when the condition breaks down, zeroth-order and first-order discontinuities can coexist on a phase 
boundary. We show with the Polyakov-loop extended Nambu-Jona-Lasinio model that such a type of 
coexistence is realized in the imaginary chemical potential region of the QCD phase diagram. We also 
present examples of coexistence of the same-order discontinuities in the real chemical potential region. 
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Exploring the phase diagram of Quantum Chromodynamics (QCD) is one of the most important 
subjects in hadron physics. Actually, many works were done so far on this subject, and it is 
expected that there appear several interesting phases in hot and/or dense quark matter; for example, 
chiral symmetry broken and restored phases, confinement and deconfinement phases, two-flavor 
color superconducting and color-flavor locked phases, and so on; for example, see Ref. [1] and 
references therein. These phases are characterized in terms of some exact or approximate order 
parameters such as the chiral condensate, the diquark condensate, the Polyakov loop, and so on. 
Therefore, correlations among these order parameters are to be investigated. In particular, the 
relation between orders of their discontinuities is essential. It was proven by Barducci, Casalbuoni, 
Pettini and Gatto (BCPG) [2] that different first-order phase transitions take place simultaneously. 
The theorem corresponds to a generalization of the Clausius-Clapeyron relation. 

Studying these correlations directly in QCD is desired, however, in the finite chemical poten- 
tial region, lattice QCD is still far from perfection because of the sign problem; for example, see 
Ref. [3] and references therein. So the phase diagram was investigated with effective models. Re- 
cently, importantprogress was made by the Polyakov 
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23(1 . This model can de- 



scribe the chiral, the color superconducting and the confinement/deconfinement phase transitions. 

Figure Q] shows the phase diagram in the chiral limit predicted by the two-flavor PNJL model 
without diquark condensate; the details of the calculation will be shown latter. The diagram is 
drawn in the /i 2 -T plane, where T stands for temperature and \i for quark chemical potential. 
The solid and dotted curves represent first- and second-order chiral phase transitions, respectively. 
In this paper, when an order parameter has a discontinuity in its value (zero th- order), we call it 
the first-order phase transition. Meanwhile, when an order parameter has a discontinuity in its 
derivative (first order) and its susceptibility is divergent, we refer to it as the second-order phase 
transition. However, our discussion is mainly concentrated on the relationship between zeroth- 
and first-order discontinuities. 

On the solid curve between points C and D, two zeroth-order discontinuities emerge simulta- 
neously in the Polyakov loop and the chiral condensate. This is a typical example of the BCPG 
theorem. As an interesting fact, on the dotted curves, two first-order discontinuities take place 
simultaneously in the chiral condensate and the Polyakov loop. This implies the BCPG theorem 
on the zeroth-order discontinuity of order parameter can be extended to the case of the first-order 
one. As another interesting point, on the dashed curve between points A and B, a first-order 
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discontinuity of the chiral condensate coexists with zeroth-order discontinuities of quark number 
density and other 9-odd quantities, where 9 = —i/i/T. On the dashed curve moving up from 
point B, furthermore, the quark number density still has a zeroth-order discontinuity, although the 
chiral condensate is always zero. Thus, the relation between orders of discontinuities of various 
quantities is much richer than the BCPG theorem predicts. 

In the left half plane of Fig. [fl \i is imaginary. However, the phase diagram in the region is also 
important, since in the region lattice QCD has no sign problem and then its results are available. 
Hence, the validity of the PNJL model can be tested there by comparing the model results with 
the lattice ones. Actually, it has been shown for the case of finite quark mass that the results of 
the PNJL model are consistent with those of the lattice simulations [220. Furthermore, the real fi 
system can be regarded as an image of the imaginary /x one, since the canonical partition function 
of real /x is the Fourier transform of the grand canonical partition function of imaginary /x n26A . 




Fig. 1: Phase diagram in the zx 2 -T plane predicted by the PNJL model in the chiral limit. 



The aim of this paper is to extend the BCPG theorem on the zeroth-order discontinuity of order 
parameter to the case of the first-order discontinuity, that is, we show that once a discontinuity of 
either zeroth-order or first-order takes place in an order parameter such as the chiral condensate, 
discontinuities of the same order appear in other order parameters such as the Polyakov loop. The 
original coexistence theorem of BCPG on the zeroth-order discontinuity and the present coexis- 
tence theorem on the first-order discontinuity are preserved, when the phase boundary is shifted 
in both the T and /x directions by varying values of external parameters such as the current quark 
mass; the condition will be shown later in (fT5l) and (fT6l) . In other words, discontinuities in mutually 



different orders can coexist only when the condition breaks. Such a situation is not just a trivial 
exception but a physical relevance. Actually, we will show that the situation is realized in the 



the Roberge and Weiss (RW) phase transition Ii26ll appearing in the imaginary chemical potential 



region of the QCD phase diagram, and prove from the viewpoint of the coexistence theorem that 
the RW phase transition is a family of zeroth- and first-order discontinuities. This resolution of 
the RW phase transition is a principal subject of the present paper. We present some examples of 
the coexistence by using the PNJL model for both real and imaginary chemical potential regions 
in the phase diagram. 

We begin with the grand canonical partition function 

Z(T,ii)=Trexp[-P(H-iiN)) (1) 

with a Hamiltonian of the form 

a 

where H determines the intrinsic system, A Q are external parameters conjugate to the hermitian 
operators O a and (3 = 1/T, /x is the chemical potential and N is the particle number. The thermo- 
dynamical potential Q(T, fx) is given by 

f2(T,fi) = ~\nZ(T,ii) (3) 

with V the three-dimensional volume, and the entropy density s and the particle number density n 
are also by 

where the subscript x means that x is fixed in the partial differentiation. The expectation value of 
the operator O a per volume 

o a = ^ = ^Tr{d a exp[-P(H - fiN)]} (5) 



is given by 



\OA a / T>,A' 



where the subscript A' shows that all the A except A Q are fixed in the partial differentiation. The 
subscripts of the partial differentiation will be suppressed for simplicity, unless any confusion 
arises. 
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First we recapitulate the original BCPG theorem yQ on the zeroth-order discontinuity (the 
first-order phase transition) in order to know what is assumed in the proof. The proof is made 
as follows. We start with the assumption that there appears a discontinuity in o 7 , and show that 
the discontinuity propagates to other order parameters o a t^ 7 . Hereafter, a' stands for a except 7. 
Thus, no assumption is made beforehand on the property of discontinuities appearing in o a /. The 
first-order phase transition appearing in o 7 is drawn by the solid curve (/x c ,T c ) schematically in 
Fig. |2l its typical example is the chiral transition at low temperature shown in Fig. [TJ The phase 
boundary (curve A) is shifted to curve B by taking different sets of external parameters, {A Q }#- 
The thermodynamical potentials Qi of phases i = 1 and 2 on curve A satisfy the Gibbs condition 

= f2 2 (T c ({X Q }) lf i c ({X a }),{X Q }). (7) 



Curve B 


Symmetry restored phase 


Curve A \. a t 








Symmetry broken phase 





Chemical potential 



Fig. 2: External parameter dependence of the phase boundary. Phase boundaries are projected on the \x-T 
plane from the X a -(j,-T space. 



Differentiating the thermodynamical potentials with respect to A 7 on the curve leads to 



an 1 



dT 



dT c dOx 
— - H 

:<9A,y 



d/i c 



dX, 



dT 



dT c dQ 2 



c dXj dfx 



*,dX~ 



(8) 
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where the subscript | c denotes that the quantities are evaluated at (/x c , T c ). Hence we obtain 



dT dp, 
os + -r—on. 



OX 



dX, 



(9) 



v 7 iv/vy 

where 5x = x 1 — x 2 (x = o 7 , s, n) is evaluated on the phase boundary. In Fig. [2l the correspon- 
dence between each individual point on curve A and that on curve B is not unique. This means 
that one can define an infinitesimal variation of T c in the T direction with fixed /i c , 

m 



T C (A 7 + AX 7 ) - T C (A~ 



dX~ 



AX 
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(10) 



and an infinitesimal variation fi c in the /i direction with fixed T c , 

dfi, 



ji c (X 1 + Z\A 7 ) - /i c (A 7 ) 



AX 



1 1 



(11) 



where the subscript \ x denotes that x is fixed; these variations are illustrated by the arrows in Fig. 
[2l Using these variations, one can see from © that 



5o^ 



5s r- 

Mc OA 



5n. 



(12) 



We find from 5o 7 ^ that 9T c /9A 7 | /ic and <9/i c /<9A 7 |T c are nonzero, since 5s and 5n never diverge. 

When all the Aq,'s are fixed at zero, yu c can be regarded as a function of T c : fx c = /i c (T c ). 
Differentiating © with respect to T c , one can get 



5n 



(13) 



Equation (fT3l) is the Clausius-Clapeyron relation between 5 s and 5n [27], and Eq. (fT2|) is a gener- 
alization of the relation to the case of nonzero X a . 
A relation similar to (PT21) is obtainable for a 1 : 



dT c 
oX n > 



d/j, c 



5s <n 

fi c (JA a > 



5n. 



(14) 



Here it should be noted that the curve (/i c , T c ) is defined by a discontinuity appearing in o 7 . The 
discontinuity 5o 7 ^ induces a new discontinuity 5o a / ^ through 5s ^ 0, when 



dT r 



dX n > 



^0. 



(15) 



/'c 



Similarly, the discontinuity 5o 7 7^ induces 5cv 7^ through <5n 7^ 0, when 



dX n i 



^0. 



(16) 



Thus, when the conditions (fl~5T) and (fT6l) are satisfied, two first-order phase transitions take place 
simultaneously. In other words, the discontinuity of o 7 propagates to other physical quantity o a > 
through those of s and n. The conditions mean that the phase boundary is shifted in both the T 
and \i directions in the fi-T plane by varying X a i . 

An early application of the BCPG theorem was to the case of a 2+1 flavor model in which two 



chiral condensations exist H28H . A similar situation is expected when an isospin chemical potential 
is introduced in 2 flavor models if a flavor mixing interaction is included [29]. 

Here we show an example of the simultaneous occurrence of zeroth-order discontinuities of 
order parameters by using the PNJL model in the chiral limit. The formulation and the parameter 
set of the PNJL model are given in Refs. [18, where the pure gauge part is obtained by 



reproducing lattice QCD data in the pure gauge theory [31 



32D as shown in Ref. (lfl]. In the 



present paper, we put m = with keeping other parameters unchanged. 

In the chiral limit, the chiral condensate a = (qq) is an exact order parameter of the spon- 
taneous chiral symmetry breaking: namely, o 7 = o and A 7 = m for the quark field q and the 
current quark mass m . The Polyakov loop <P is an exact order parameter of the spontaneous Z 3 
symmetry breaking in the pure gauge theory, but the symmetry is not exact anymore in the sys- 
tem with dynamical quarks. However, <P still seems to be a good indicator of the deconfinement 
phase transition. We then regard $ as an approximate order parameter of the deconfinement phase 
transition. 

Figure [3ta) represents T dependence of a, <P and the charge-conjugated Polyakov loop 3 at 
\i = 280 MeV. These are discontinuous at the same temperature T = 154 MeV. This behavior 
is consistent with the BCPG theorem that guarantees the simultaneous occurrence of zeroth-order 
discontinuities of order parameters. Here one can also find that ll dependence of <& is similar to 
that of <P. This is true for other real ji. 



The BCPG theorem does not necessarily mean that there can not exist a quarkyonic phase H33H 
defined in the limit of large number of colors as a phase that has a finite value of the baryon 
number density n but is confined. This is understandable as follows. The quarkyonic phase was 



recently investigated with the PNJL model [7, 



2li 



350 and the strong coupling QCD [|34ll . The 



PNJL analysis of Ref. 112 lh shows the simultaneous occurrence of zeroth-order discontinuities of 
a, <1> and n. However, the discontinuity in $ is only a jump from a small value to another small one, 
while that in n is a jump from almost to a value larger than the nuclear saturation density. As 
mentioned above, <P is only an approximate order parameter of the deconfinement phase transition, 
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and then such a small jump is possible. Such a small jump could be just a propagation of the 
discontinuity in a. Thus, we can not say necessarily from the small jump of <P that a first-order 
deconfinement phase takes place together with the chiral phase transition and the phase transition 
of n. A plausible definition of a critical temperature of the deconfinement phase transition is a 
temperature that gives $ = 0.5. In this definition, the deconfinement transition is crossover, and 
then the BCPG theorem is not applicable anymore. 




T [GeV] T [GeV] 



Fig. 3: The temperature dependence of the chiral condensate a, the Polyakov loop <P and its conjugate <P at 
(a) fj, = 280 MeV and (b) fi = 50 MeV in the chiral limit. The chiral condensate is normalized by the value 
at T = \i = 0. The inset figure in (b) represents <P and 3 near T = 260 MeV that is the critical temperature. 

Next, we proceed to the case that an order parameter has a first-order discontinuity. In this 
case, since a first-order discontinuity becomes continuous by a change of external parameters, the 
boundary must be defined in terms of a susceptibility as follows. Here we take the chiral transition 
in 2 flavor systems at high temperature, shown by the dotted curve in Fig.[H as a typical example: 
namely, A 7 = m and o 7 = a = (qq). The second-order chiral phase transition at m = 
becomes crossover whenever m is finite B22I1 . It is possible to define the phase boundary of such 
a crossover with the chiral susceptibility \<j — —da/dm so that the T dependence of % becomes 
maximum on the boundary. This definition works also in the chiral limit, although the maximum 
is infinity. In this case, curves A and B in Fig. [2] are reinterpreted as phase boundaries so defined 
and correspond to the chiral limit (m = 0) and to the case of small m , respectively. Curve 
A can move continuously and reach curve B by varying m from to a finite value. The chiral 
susceptibility Xa is divergent on the boundary A, since the chiral phase transition is of second 
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order there. 

Now, we consider curve A defined above and its vicinity. The system concerned has no zeroth- 
order discontinuity, 5s = 5n = 0. Differentiating 5s = with respect to A 7 on the boundary 
(/i c , T c ) leads to 

Using the relation ds/dX 1 = —do a /dT and the variations in the T and /i directions mentioned 
above, one can get 

M c \<rr) ox. 



dT J <9A 7 



S(¥). (18) 



t c \dfiJ' 



Taking the same procedure for 5n = 0, one also obtains 

\ d[i) <9A 7 nc^xdr) 8X, f T c \d/j)' 

Other order parameters o a / satisfy the same equations as (TTBl and (fT9l ). Note that all the equations 
are evaluated in the chiral limit A 7 = m = 0. It is found from (TTSl and (fT9b for o 7 and the 
corresponding equations for o a > that discontinuities 5(do 7 /dT) ^ and 5(do 7 /dfi) ^ induce 
new ones 5(do a > /dT) ^ and 5(do a //dfi) ^ 0, when the conditions (TT5T) and (fT6l) are satisfied. 
Thus, two first-order discontinuities of order parameters can coexist under the conditions (fT5l) and 
(TT6l) . Furthermore, it is found from (fTBl and (fT9b that 5{do 1 /dij) is not zero whenever 5{do 1 /dT) 
is not zero, because of dn/dT = ds/dfi. Accordingly the first-order discontinuity of o 7 emerges 
in both do^/dT and do^/dfx. 

Here we show an example of the simultaneous occurrence of two discontinuities by the PNJL 
model in the chiral limit. Figure Ob) represents T dependence of a, <P and 3 at = 50 MeV. 
Obviously, $ and <P are not smooth at T c = 260 MeV. In the inset figure of Fig. Ob), the solid 
curves show <P and <P near T c = 260 MeV, and two dotted lines do tangential lines of the solid 
curves at T = T c — 0. The deviations between the solid curves and the corresponding dotted lines 
indicate that <P and § are not smooth at T c = 260 MeV. Thus, da/dT, d<P/dT and d$/dT are 
discontinuous at the same temperature, as expected from the coexistence theorem on the first-order 
discontinuity of order parameter. 

As shown in Fig. Ob), 5(do/dT) diverges at T = T c , because da/dT\ T =T c -o = oo and 
da / dT\ T=Tc+0 = 0; note that a < 0. Similar divergence is also seen on the second-order chiral 
phase transition line (dotted curves) in Fig.Q] This divergence indicates from (TTBl that dT c / dmo]^ 



and/or 5(ds/dT) diverges on the second-order chiral phase transition curve. If S(ds/dT) is infi- 
nite there, the divergence will propagate to other quantities 8{do a > / dT) when the condition (fl~5l) is 
satisfied. As shown below, this is not the case of the second-order chiral phase transition. Figure 0] 
presents T dependence of ds/dT, the chiral susceptibility Xa and the Polyakov-loop susceptibility 
at 11 = m = in panel (a) and m dependence of T c at /i = in panel (b); definitions of Xa 
and x$& are shown below. As shown in panel (a), x<j ( tne dotted curve) diverges at T = 261 .4 MeV, 
but ds/dT (the solid curve) has a finite gap there. Meanwhile, panel (b) shows that the gradient 
dT c /dm is divergent at m = 0. In the present case, thus, the divergence in 8(da/dT) does not 
propagate to other quantities 5{do a i / dT). Accordingly, the coexistence of first-order discontinu- 
ities of order parameters takes place, but the coexistence of second-order phase transitions does 
not occur necessarily, because there is a possibility that 9T C /9A 7 | /Jc and <9/x c /<9A 7 |t c diverge. In 
other words, when dT c / '<9A 7 |^ c and dT c /d\ a i are nonzero and finite, or when <9/x c /<9A 7 |t c and 
dfx c /d\ a '\T c are nonzero and finite, the coexistence of second-order phase transitions takes place. 
Susceptibilities Xij °f °"> ^ and § can be written as Jol[l5, 



18] 



(20) 



where 



K 



/ 



V 



j3 d 2 Q 
4G%A da 2 

13 d 2 Q 
' 2G B A 2 d<Pda 

(3 d 2 Q 



13 d 2 Q 
2G S A 2 dad$ 

13 d 2 Q 
A3 d<P 2 

13 d 2 Q 



13 d 2 Q \ 
2G a A 2 dad<P 

(3 d 2 Q 
A 3 d<Pd<P 

(3 (Pn 

A% d<P 2 ) 



(21) 



' 2G S A 2 d<Pdcr A^ d<Pd<P 

is a symmetric matrix of curvatures of Q and (K' 1 )^ is an element of the inverse matrix 
K~ l . In the chiral limit, Q is invariant under the transformation a — > —a \W\ and hence a-even. 
In the case that the chiral phase transition is the second order, as shown in Fig. [3l Q becomes 
minimum at a = when T > T c . Therefore, K a <p and K c § are zero at T > T c , because they are 
cr-odd. In this situation, the Xij at T > T c are reduced to 

1 



(22) 



where 



13 d 2 n /3 d 2 a 



A 3 d<P 2 A 3 d<Pd<P 
f3 d 2 (2 13 d 2 Q 



(23) 



A 3 d<Pd<P A 3 d<P 2 

Thus, the susceptibilities of <& and $ are decoupled from that of a. In particular at T = T c , the 
curvature K aa is zero and then Xa is divergent, while at T > T c the curvature K a(T is positive 
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and then Xa is a positive finite value. This divergence makes no influence on other susceptibilities 
Xij = $■,$), since K a $ and K a § are zero; see Refs. 1124 12511 for the detail. Figure Ufa) is a 
typical example of this situation; Xc has a divergent peak, while x$$ does not. 
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Fig. 4: Properties of thermal quantities at fi = 0: (a) T dependence of chiral and Polyakov-loop suscepti- 
bilities and ds/dT at = and (b) quark-mass dependence of the critical temperature. In panel (a), the 
solid, dashed and dotted curves represent ds/dT, and Xa, respectively. The chiral susceptibility x<j 
has a divergent peak at T = T c = 261.4 MeV. 

In the imaginary chemical potential region (jj 2 < 0), as shown by the dashed curve between 
points A and B of Fig. [fl there coexists a zeroth-order discontinuity of quark number density n 
and a first-order discontinuity of chiral condensate a. The coexistence is consistent with the proofs 
mentioned above, as shown below. This is the principal subject of the present paper. 

It is convenient to introduce a new variable 9 = —ifx/T instead of ji. The conditions (TTBI) and 
(TT6l) are then changed into 

dT r 



d\ a i 



d\ n > 



+ 0, (24) 

+ 0. (25) 

As shown later in Fig. [71 the coexistence of 5n ^ 0, 5a = and 5(da/d9) ^ always appears 
on vertical lines 9 = (2k + l)n/3 in the 9-T plane, where k is an integer. This indicates that 
d9 c /dm = and then the condition (1251) breaks down. Taking the new variable 9 also changes 
© into 

r~ &T a d0 o 

oa = - — os + - — on = - — os (26) 

om orriQ om 
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with 



tin 
ar 



n 



9s )t,x 



inT, 



a 



\drm) 6TX , 



(J. 



(27) 
(28) 



where use has been made of d9 c /dm = in the second equality of (l26t . Thus, even if 5n is not 
zero, one can keep 5a = in (l26l) when 5s is zero. 



The discontinuity 5h ^ was first pointed out by Roberge and Weiss (RW) [26], and of- 
ten called the RW transition. Here, we consider how the discontinuity 5h ^ influences 
other order parameters o a . In this case, curve A in Fig. [2] is defined by the discontinuity 
5h 7^ 0. The quantity 5h is a function of T c , 9 C and X a , but 9 C does not depend on X a ; namely, 
5h = —f(T c ({X a }),9 c ,{X a }). Differentiating 5h + / = with respect to X a leads to 



5 



dn 



+ 



dT c 



5(t) + 



df 



+ 



dT 



c d\ n 







(29) 



because of d9 c /d\ a = 0, where the subscript | £ denotes that the quantities are evaluated at 
Using dn/d\ a = —do a /d9 and taking the variation in the 9 direction with fixed T c , one can obtain 

■do a \ df 



5 



89 I d\ a 
Taking the same procedure for 5s = leads to 

' do n \ 89, 



(30) 



5 







(31) 



dT J dXc_ 

because of d9 c /dX a \T c = 0. Equations (|30l ) and (|3TT) indicate that order parameters are discontinu- 
ous in do a /89 but not in do a / dT. This property is different from that of the ordinary second-order 
chiral phase transition, shown by the dotted curve of Fig. [0 that is discontinuous in both do a /d9 
and do a jdT . The coexistence of the zeroth-order discontinuity of n and the first-order one of a 
is originated from the fact that the RW transition line is vertical in the 9-T plane and does not 
move in the 9 direction by changing the external parameter X a = m . The influence of the RW 
discontinuity of n to the Polyakov loop is discussed in the following. 

In the imaginary /i region, physical quantities have a periodicity of 27r/3 in 9, when these are 
invariant under the extended Z 3 transformation B22I1 . 



JziO 



$(9) -> $(0)e 



$(9) -»■ $(6)e 



3 



(32) 
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This is called the RW periodicity Q26Q . The thermodynamical potential i?pNjL and a, s and n 
are invariant under the extended Z 3 symmetry, but <P and <P are not [12211 . However, this can be 
cured by introducing the modified Polyakov loop \P = <P exp(i9) invariant under the extended Z 3 
transformation. We then consider a period < 9 < 2n/3 without loss of generality. 

Figure [5] shows 9 dependence of the chiral condensate a and the imaginary part of n, Im[n], at 
T = 300 MeV; note that n is pure imaginary for imaginary /i by definition. The chiral condensate 
has a cusp at 9 = n/3, while n is discontinuous there. Thus, the first-order discontinuity of a and 
the zeroth-order discontinuity of n coexist, as predicted above. 



0.02r 



1 0.01 

D 
I 



0.05 



(a) 



CD 



2 

e/(ji/3) 



-0.0! 



(b) 




2 

6/(7i/3) 



Fig. 5: The 9 dependence of the chiral condensate and the imaginary part of the quark number density at 
T = 300 MeV in the chiral limit. 



Next, we consider the relation between the discontinuity of n and the Polyakov loop transition 
by using Re[^] = {${9) + W(9))/2. In the PNJL model, &(9) and &(0) are treated as classical 
variables, and it is found from the expression for i7 PNJL in (13) of Ref. fl 1 8f1 that &(0) is the 
complex conjugate of #(9) for the case of imaginary //. Figure |6£ a) shows 9 dependence of the 
real part Re[^] at T = 300 MeV. There appears a first-order discontinuity also in Ref!^] on the line 
9 = 7r/3, as expected from (|30l . 

Finally, we consider the imaginary part of Imf^] = (\P(9) — &(6))/2i. This is also real, but 
#-odd (odd under the interchange of 9 <-> — 9 ), because &(9) = ^(—9) j22\\ . One can not use 



A a Im[?^] as a source term X a O a , since it breaks ^-evenness, i?pNjL(^) = ^pnjl(— 9), that is the 
charge-conjugation symmetry of ^pnjl [36]. To avoid this problem, we introduce a source term, 
A a sin(3^)Im[<^], designed to keep ^-evenness and the RW periodicity. This is just an example of 
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operators having the two properties. For this source term, (|3Q|) is reduced to 

1 df 



3d\ r 



(33) 



Thus, 5(Im[<? r ]) is finite on the RW phase transition line 9 = ir/3 because of df /d\ a \c ^ 0. This 
indicates that the zeroth-order discontinuity of n induces that of Imflf] as shown in Figlfjb). 
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Fig. 6: The 6 dependence of !f at T = 300 MeV in the chiral limit:(a) the real part and (b) the imaginary 
part. 

Throughout all the analyses, we can conclude that the zeroth-order discontinuity of a 0-odd 
quantity n induces zeroth-order ones in 0-odd quantities and simultaneously does first-order ones 
in 0-even quantities; see Ref. [22] for the proof of the even/odd property of n, a, Re [if'], Im[?f], |!f| 
and arg[?f]. As shown in Fig. [5] da/ 09 is finite on both the sides of the critical chemical potential 
9 C . This means that the chiral susceptibility x is finite. Hence, there is no second-order phase 
transition on the RW line. Therefore, the RW phase transition is a first-order phase transition and 
a family of zeroth- and first-order discontinuities. 

Figure [7] shows the phase diagram on the 9-T plane that corresponds to the /j 2 < part of Fig. 
[Q On the dashed line between points A and B, the RW phase transition mentioned above emerges. 
The transition comes out also on the dashed line going up from point B, although a is zero there 
and then no discontinuity takes place in a. The dotted curves represents ordinary chiral phase 
transitions of second order. 

To summarize, we showed that once a zeroth- or first-order discontinuity takes place in a quan- 
tity o 7 , discontinuities of the same order emerge in other quantities o Q ^ 7 , if the conditions (fT5l) and 
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Fig. 7: Phase diagram on the 6-T plane predicted by the PNJL model in the chiral limit. 

(fT6l) are satisfied, that is, if the phase boundary is shifted in both the directions of T and \i in the 
T-fj, plane by varying values of external parameters X a conjugate to o a . This coexistence theorem 
is an extension of the BCPG theorem on the zeroth-order discontinuity of order parameter (the 
first-order phase transition). When the conditions break, first- and second-order discontinuities 
can coexist on the same phase boundary. The RW phase transition in the 6-T plane, composed of 
zeroth-order discontinuities of 0-odd quantities and first-order discontinuities of ^-even ones, is a 
typical example of the coexistence of zeroth- and first-order discontinuities. The RW phase tran- 
sition line is vertical and does not move in the 6 direction, even if any external parameter varies. 
Thus, the shape of the phase boundary and its variation with external parameters are essential in 
determining which type of coexistence is realized. 

This work has been supported in part by the Grants-in-Aid for Scientific Research (18540280) 
of Education, Science, Sports, and Culture of Japan. 
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